Quantum insulating states of F = 2 cold atoms in optical lattices 
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In this Letter we study various spin correlated insulating states of F = 2 cold atoms in optical 
lattices. We find that the effective spin exchange interaction due to virtual hopping contains an 
octopole coupling between two neighboring lattice sites. Depending on scattering lengths and num- 
bers of particles per site the ground states are either rotationally invariant dimer or trimer Mott 
insulators or insulating states with various spin orders. Three spin ordered insulating phases are 
ferromagnetic, cyclic and nematic Mott insulators. We estimate the phase boundaries for states 
with different numbers of atoms per lattice site. 



Recently, the scattering lengths of Kb atoms in the F=2 
manifold in optical lattices have been studied experimen- 
tally These results, together with earlier theoretical 
estimates of scattering lengths offer critical informa- 
tion about how atoms interact in the F — 2 manifold 
and are vital for the further understanding of their mag- 
netic phases. Investigation of cold atoms with high spins 
in optical lattices not only improves our understanding 
of fundamental ideas in quantum magnetism, but might 
also lead to major breakthroughs in quantum informa- 
tion storing and processing 0, [j] . A few impressive theo- 
retical efforts were already made to understand BECs of 
F = 2 sodium and rubidium atoms in optical traps @, Q • 
Depending on the two-particle scattering lengths in total 
spin F = 0, 2, 4 channels, the ground states of cold atoms 
can be either polar, cyclic or ferromagnetic condensates. 
The existence of various spin correlated states of cold 
atoms in optical lattices has also been suggested @, [3 ■ 

In this Letter we shall address the following question: 
What are the possible insulating states of F = 2 cold 
atoms in optical lattices? We first summarize the main 
results of our investigation, a) When optical potentials 
are strong and exchange interactions are asymptotically 
weak, for certain numbers of particles per lattice site 
the ground states are rotationally invariant Mott insu- 
lators of dimers (i.e., singlets of two F = 2 atoms) or 
trimers (i.e., singlets of three F = 2 atoms), b) Decreas- 
ing the optical potential depth, depending on scattering 
lengths, results in various spin ordered insulating states 
such as cyclic, nematic and ferromagnetic insulators, c) 
The boundaries between different insulating phases are 
estimated (See Fig.([T])). 

To obtain these results, we consider F — 2 cold atoms 
which have two-body scattering lengths, (Zf=o,2,4j m 
three channels with total hyperfine spins F = 0, 2, 4. The 
two-body contact interaction is X)f=o i 4 9f~Pf where 
gp = 4:irh 2 ap/M and Vf are the projection operator 
which projects out total spin F states of two atoms. 
Alternatively this can be rewritten in terms of a + 
b/2{F 2 - 12) + 5cP where F is the total spin oper- 
ator. Three interaction constants were calculated pre- 
viously, a = (4cj 2 + 3g 4 )/7, b = -(g 2 - 9i)/7 and 



c = (.90 - 54)/5 - 2{g 2 - g 4 )/7 @, 0. To study the 
states in optical lattices, we find it is convenient to intro- 
duce a tensor representation 4b> a, f3 = x,y, z for F = 2 
states: 
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where V4if?> rap — 0, ±1,±2 are the usual creation op- 
erators for F = 2 particles. The tensor operators 4 a p 



are symmetric and traceless, that is ijy 



4 0a 



and 



trij}^ — 0. We find this a convenient labelling of the 
five F — 2 states. The commutation relations are 



bl>af3,4 a >{3>] = Saa'Spp' +5 a 0'5f3 a > - -S a ffS a 'ffi. (2) 

The construction of rotationally invariant operators in 
this representation is straightforward. For instance, the 
number operator p, the dimer or singlet pair creation 
operator 2?t and the trimer or singlet of three atoms cre- 
ation operator T' are 
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(3) 

where ip^ represents the tensor. The total spin operator 

F a , a — x, y, z is 

F a = -iea^^Wr) [F*iFp\ = (4) 

Finally, F a commutes with rotationally invariant op- 
erators such as X?t and T^; generally one finds that 
[F a , tr(ip') n ] = 0. For F — 2 cold atoms in optical lat- 
tices, we employ the following Hamiltonian 
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+ 5CL V l V * ~ 1 E ^taP^Pa + k - C -)- ( 5 ) 
fe <fci> 

Here fc is a lattice site index and < kl > are neighboring 
sites, t is the one-particle hopping amplitude, cll^Li^l 
are, respectively, proportional to a, 6, c introduced earlier; 
glC&l, cl) = ^(6, c) / dr^Qtpo) 2 and ^ is the on-site lo- 
cal ground state wavefunction. plays the role of 
projection operator Vq discussed before. For most cold 
atoms, the absolute values of &l, cl are much smaller 
than ol and we work in this limit unless specified. When 
the optical potential depth is small and t is much larger 
than ol, cold atoms effectively are weakly interacting and 
condense. The results obtained using Eq.© are identi- 
cal to the ones previously obtained and we do not present 
them in this Letter. Here we focus on the physics of in- 
sulating states when t is much smaller than a^. 

The effective low energy Hamiltonian for a Mott insu- 
lator with M(> 1) atoms per lattice site is 

« = yE(^ 2 - 6 ^) + 5c ^ E v l v >< 

k k 

-Jo E (9k;aP,a>p>QhaP,a'P +h.c)\ 
<kl> 

f 2 

— trtl)lipk(6 a i a 6ppi + 5 a p>5f3 a > - -S a p5 a >p/) (6) 

where Jo = t 2 /cl- The virtual hopping between two 
neighboring sites results in an exchange interaction of 
the form of octopole coupling as shown above. This is 
valid for all Mott states with the number of particles M 
larger than one. An exchange interaction of quadrupole 
type has been derived recently for F = 1 cold atoms in 
optical lattices @, H, E3] • 

The spin-ordered Mott insulating states under consid- 
eration (with the crystal translational symmetry) are 

i*>=n*^i° =5; (?) 

X is a traceless symmetric complex tensor and M is the 
number of atoms per site. We calculate the x _ dependent 
part of the mean field energy per particle, 

16 

Emi = M - 1 )A X ,X* - —J ex M; 

5 

A(x, X*) = 4b L Tr[ Xl x*] 2 + 2c L tr X 2 Tr X * 2 (8) 

Here J ex — zJo, z is the number of neighboring lattice 
sites. Minimization of this energy subject to the con- 
straint of the normalization of x m Eq. Q) can be carried 
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FIG. 1: (Color online) Schematic of phase diagrams for spin 
correlated Mott states of F = 2 cold atoms in optical lat- 
tices. In (a)-(d), the phase diagrams are plotted in a plane 
of two interaction constants bL, cl for different M (the num- 
ber of particle per lattice site). Only insulating states are 
shown. a_t, are defined in terms of scattering lengths after 
Eq. ([5]). In (e) and (f), the vertical axis is the hopping t and 
the horizontal axis is the chemical potential ^t; the numbers 
along the ^i-axis indicate the number of particles per lattice 
site in Mott states (states inside lobes). Outside the lobes 
are superfluids which we do not discuss in this article, (a) for 
M = 3k = 2k'; (b) for M = 3k # 2k' (c) for M = 2k' 3k; 
(d) for M ^ 3k ^ 2k' or for M = 1. (e) for -c L > b L > 0; 
(f) for cl S> bL > 0. N,C,F, D, T stand, respectively, for 
nematic, cyclic, ferromagnetic, dimer and trimer Mott insu- 
lating states; In (b) and (c), there are no phase transitions 
across dashed lines (guide for eyes only). Phase boundaries 
(thick blue lines) are estimated by applying Eq. ([8) and com- 
paring the energies obtained in Eq. (|12lll8p . 

out using the method of the Lagrangian multiplier. We 
provide the results here: 

1) bL > and cl > 0. The traceless tensor x satisfies 

tr X 2 = \x,X*] = 0- (9) 

There are two discrete root solutions XR± which have 
only diagonal components, Xxx = l/v 7 ^, Xyy = 
exp(±i27r/3)/v / 6 and Xzz — exp(±£47r/3)/-\/6 and Xap = 
when a [3. They are also invariant under a cyclic 
subgroup C3 of the SO (3) rotation group. The full set 
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of solutions can be obtained by applying SO (3) rotation, 
and U(l) gauge transformation. These solutions corre- 
spond to cyclic insulating states. 

2) &l > Cl/4 and cl < 0. Up to a phase, \ satisfies 



tr X ' 



1 

2 ,X 



(10) 



For diagonal matrices which satisfy the above conditions, 
the matrix elements have to fall on an elliptical curve; 
that is, xL+ XxxXyy +Xy V = V 4 , tr\ = and Xa/3 = 
if a ^ /3. The full set of solutions can be obtained by ap- 
plying SO(3) and U(l) rotation. These states are char- 
acterized by traceless real symmetric tensors (up to a 
phase) and < F a >= 0, a — x, y, z; moreover, the corre- 
sponding expectation value of the ncmatic order operator 
ipa^ipjp — l/35 a ptrip*ip is nonzero. They represent ne- 
matic insulating states (see more discussions in RefE3). 
3) &l < and < c/,/4. x satisfies 

tr X 2 = 0, X*XX* ~ X*X*X - tr\x, X*?X* = (11) 

and [x, X*] 0- An example of solutions is Xza = 0(a = 
x,y,z) and Xxx = ~Xyy = -iXxy = 1/V8- This corre- 
sponds to a ferromagnetic insulating phase. 

So three spin ordered Mott insulating states are fer- 
romagnetic, nematic and cyclic states (indicated by sub- 
scripts F, N, C respectively) ; the corresponding energies 
per atom are 



E F = 2(M-l)&i 
E N = (M-1)^-- 
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When 6l and cl are both positive and the exchange 
coupling is zero, the ground state of the Hamiltonian is 
still an M-atom spin singlet but has to be annihilated by 
the dimer operator T>. For M = 3, such a state can be 
created by the trimer creation operator T; the wavefunc- 
tion and its energy are 



l T >=n^M^) 3 |0>,£r 



-3b, 



(15) 



MJ ex {\2) 



Note that V\T >= and < T\V^V\T >= 0. 

For more than three particles, one can show that sin- 
glet states of the form trip n with n > 3 can always be 
rewritten in terms of two fundamental singlet creation 
operators: trip 2 for a dimer and trip 3 for a trimer. Gen- 
erally, singlet states for M = 2ki + 3fc 2 atoms can be 
expressed in terms of (trip'ip^) kl (trip'ip^ip') k2 . For in- 
stance, direct calculations show that tr(ip') 4 \0 > is iden- 
tical to (trip^ ip^) 2 \Q > and ir(-0 t ) 5 |O > is identical to 
tr(ip') 2 tr(ip^) 3 \Q >. To construct an M-particle singlet 
without singlet pairs, the only candidate to consider is 
(trip'ip^ip^) k \0 > since all other states explicitly involve 
(trip'ip') n , the dimer creation operator. Therefore, for 
an arbitrary M, a 3k-atom singlet without dimers can 
be written as \T >= H k V T (k)(tr(ipi) 3 ) k ; ( up to the 
normalization factor) here Vr(k) is a projection oper- 
ator which further projects out states without dimers 
at site k, i.e. VkPrik) = 0. Taking into account that 
F%\T >— T>k\T >— 0, we find that the energy per atom 
for a trimer state is independent of M, Et = — 3&j> It 
has a lower energy than the cyclic state when 



In addition, one should also take into account rotationally 
invariant insulating states. The physics of these states 
depends on the numbers of particles per lattice site and 
below we discuss different situations separately. 

i)M = 2k = 3k' (k, k' are integers). When hh and — cl 
are positive and J ex approaches zero, the Hamiltonian 
has the lowest energy eigenstate which is rotationally in- 
variant and has a maximal number of singlet pairs, or 
dimers. It has the following wavefunction 



\D>-- 
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(^ + 3)y-36 L (13) 



The state \D > (up to a normalization factor) is an in- 
sulator of singlet pairs or dimers; following Eq.(|12p and 
Eq. (Ti~3")) . it has a lower energy than the nematic state 
when the exchange interaction is small, i.e., 

Jex < J D>J C D = J^ ( 3& i - 2C L) ■ (W) 



In 

< 



ex u ' u 16M V " ^ 

dimer Mott states, all atoms are paired in singlets and 
D\V^V\D >= M(M + 3)/10. 



16M 



(16) 



By further comparing Ed,t and Ec.n.f we obtain the 
mean field phase diagram; we find the mean field phase 
boundaries between rotationally invariant trimer and 
dimer states and three spin ordered phases, taking into 
account Eq. ([T2]) . Eq. ([13"]) and Eq. ([l5]l . The results are 
summarized in Fig.(fl]a). The boundary between T- and 
D-phase is cl — 0, C- and F-phase is &l = 0, T and 
C-phasc is = 16MJ ex /15, between D- and A^-phase 
is &l = 2cl/3 + 16MJ ex /15, between N- and F-phase 
is bi = Cl/4, and finally between D and F-phase is 
b L = (M + 3)cl/(4M + 2) + 16MJ ex /5(2M + 1). 

ii) M = 3k = 2k' + 1. Following the general construc- 
tion outlined in i) , the trimer states are ground states in 
optical lattices when cl and &l are both positive and J ex 
is approaching zero. However, when cl is negative, 
is positive, as J ex approaches zero, at each lattice site k' 
singlet pairs or dimers are formed leaving the last atom 
unpaired. So the system effectively is a lattice with an 
5 = 2 spin at each lattice site. The situation is similar 
to what happens to F = 1 atoms in optical lattices with 
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an odd number of atoms per site [lfj. Following the dis- 
cussions on the case for one atom per lattice (M = 1) 
below, we find that the ground state remains to be ne- 
matically ordered down to J ex = as far as cl is negative 
and bh is positive. Only four phases (trimer Mott insu- 
lators and spin ordered cyclic, ferromagnetic, nematic 
Mott states) appear in the phase diagram as shown in 
Fig.(fT]b). The phase boundary between T- and iV-phase 
is &x = (1 — M)cl/6 + 16MJ ex /15 which approaches 
cl = as M becomes infinity. The other phase bound- 
aries are the same as in Fig. (JT]b) discussed before. 

iii) M = 2k = 3k' ± 1. When cl is negative and \>l 
is positive and the exchange interactions are small, the 
ground states are dimer insulating states. When cl and 
&l are positive and J ex approaches zero, at each lattice 
site kl trimers are formed; the last atom (for M = 3fc'+l) 
or two (for M = 3k' + 2) which do not participate in 
the trimer formation are in an F — 2 spin collective 
state. Effectively at each site only the F = 2 state is 
relevant for spin ordering and the system is again equiv- 
alent to a lattice spin Model for S = 2 spins. This limit 
is equivalent to M = 1 case studied below and we find 
that the cyclic state remains to be the ground state as 
J ex becomes zero, provided that &z, and cl are positive. 
Four phases that appear in the phase diagram are dimer 
Mott insulators and spin ordered cyclic, ferromagnetic, 
nematic Mott states, as shown in Fig.fHc). The bound- 
ary between D and C is b L = l6MJ ex /15 + (M + 3)/6c L ; 
again it approaches cl — as M becomes infinity. The 
other boundaries are identical to those in Fig. (fl] a) . 

iv) Finally, we would like to address the special case 
of M = 1. The effective Hamiltonian derived before is 
only valid when M > 1. Here we directly calculate the 
exchange interaction between two adjacent F — 2 spins 
due to virtual hopping and obtain the following Hamil- 
tonian 

u <kl> F=0,2,4 

Here VfQzI) projects out a total spin F state of two 
atoms on neighboring sites < kl >. a is defined as 
a = (4d2 + 3ai)/7. For two adjacent atoms with total 
spin F = 1,3, the exchange interaction vanishes. We 
calculate the energy per particle using the wavefunction 
specified by Eq. ([7]) with M = 1 and derive an expression 
for the energy as a function of x- The x-dependence of 
the energy per atom can be shown identical to that in 
Eq.©. We again obtain three phases and the energies 
per atom of cyclic, nematic and ferromagnetic states are 

Ep a E c 3a 45, 
Jo 04 ' J 7a 4 7a2 ' 

E N _ 18a 2a a_ , 

Jo 35a4 7a2 5ao 



When the difference between scattering lengths is small, 
one can easily verify that the phase boundaries of these 
three insulating states are identical to those of three con- 
densates discussed previously. The results are summa- 
rized in Fig.([T]d). The boundaries between F and C, C 
and N and N and F are, respectively, ol = 0, Cl = and 
bL = Cl/4 ill|. Finally we present the results in terms 
of hopping t and chemical potential jj, of F = 2 atoms 
in two parameter regions: 1) — cl » 6l > ( Fig. ((It)); 
and 2) cl » &l > (Fig.((U) ); the spin ordering exhibits 
mod 2 (Fig.(fTt)) and mod 3 (FigJTf) behaviors. 

In conclusion, we find that F — 2 cold atoms in op- 
tical lattices in addition to forming spin ordered Mott 
states, can also be in states that only consist of dimers or 
trimers. These results are relevant to Kb and Na atoms 
in optical lattices. According to earlier estimates^], 
87 Kb atoms fall nearly on the cl = line (&£ > 0), a 
phase boundary between the cyclic and nematic insulat- 
ing phases. For 23 Na, bL = — 4.61c^ > 0; the insulating 
state of 23 Na atoms is nematic when J ex is large; for 
M = 2k as the optical potential decreases, a quantum 
phase transition occurs and the ground state becomes 
a dimer Mott state. However, the error bars in esti- 
mates of scattering lengths might lead to uncertainties 
in determining ground states. Finally, dimer or trimer 
Mott states have fully gapped spin excitations and are 
robust; spin ordered states in 3D optical lattices are sta- 
ble against quantum fluctuations (see similar discussions 
in Ref.[I(|). In a finite trap because of coexistence of 
Mott states with different occupation numbers flij]. the 
spin order spatially alternates between dimer (trimer) 
and nematic (trimer) ones. This work is supported by 
the office of the Dean of Science, UBC, and NSERC, 
Canada. FZ currently is an A. P. Sloan fellow. 
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